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A consistent first-order model for relativistic heat flow 

Nils Andersson & Cesar Lopez-Monsalvo 
Mathematical Sciences, University of Southampton, UK 

This paper revisits the problem of heat conduction in relativistic fluids, associated with issues concerning both 
stability and causality. It has long been known that the problem requires information involving second order 
deviations from thermal equilibrium. Basically, any consistent first-order theory needs to remain cognizant 
of its higher-order origins. We demonstrate this by carrying out the required first-order reduction of a recent 
variational model. We provide an analysis of the dynamics of the system, obtaining the conditions that must 
be satisfied in order to avoid instabilities and acausal signal propagation. The results demonstrate, beyond any 
reasonable doubt, that the model has all the features one would expect of a real physical system. In particular, we 
highlight the presence of a second sound for heat in the appropriate limit. We also make contact with previous 
work on the problem by showing how the various constraints on our system agree with previously established 
results. 



I. CONTEXT 



Relativistic thermodynamics continues to provide interesting challenges, in particular in the context of dissipative and nonlin- 

|j;^ ear phenomena. The issues involved range from direct applications in various areas of physics to fundamental problems like the 

nature of time (visavi the second law of thermodynamics) and the formation of structures at nonlinear deviations from thermal 

equilibrium. Much recent work has been motivated by the modelling of complex astrophysical systems, like neutron stars 111], 

I— I and cosmology |2]. There has also been a resurgence of interest in dissipative systems in the context of colliders like RHIC at 

O Brookhaven and the LHC at CERN OJ-qII- These latter developments, which have to a large extent been driven by the need to 

^ ' understand the dynamics of a hot quark-gluon plasma, are often linked with underlying principles like the AdS/CFT conjecture 

?H and holography [6]. Even though the problem dates back to the origins of relativity theory, it remains (in a slightly different 

, ^^ guise) at the forefront of modern thinking. 

According to the established consensus view, one must account for second-order deviations from thermal equilibrium in order 
T— I . to achieve causality and stability. This is certainly the lesson from the celebrated work of Israel and Stewart Q aL see |I9!- 
E> ' [nl] for recent work on the problem. We have recently revisited the key points in the context of heat conduction III2I1 . taking 
y^ a multi-fluid prescription based on Carter's convective vaiiational formulation for relativistic fluids |13] as our starting point. 
This is a mathematically elegant approach that has the flexibility required to account for the physics that we need to consider. 
A particularly appealing feature of the variational approach is that, once an "equation of state" for matter is provided, the 
theory provides the relation between the various currents and their conjugate momenta |1]. The variational analysis leads to 
C^ a second-order model which has the key elements required for causality and stability, in particular, it clarifies the role of the 
inertia of heat (e.g., the effective mass associated with phonons). This effect enters the model in an intuitive fashion in terms 
of entrainment between the matter and heat fl4]. As demonstrated by Priou [15] some time ago, the final variational model 
is formally equivalent to the Israel-Stewart construction. This exercise demonstrates clearly that the relaxation associated with 
causal heat transport is determined by the thermal inertia. At the end of the day, the theoretical framework becomes rather 
intuitive and the physics involved seems natural. 
^ Does this mean that no troublesome issues remain in this problem area? Not quite. First of all, it is clear that the need to intro- 

Cd duce additional parameters (e.g., the relevant relaxation times) and keep track of higher order terms (fluxes of the fluxes etcetera) 
make actual applications rather complex. Secondly, we are not much closer to considering systems that deviate significantly 
from equilibrium, such that there is no natural "small" parameter to expand in. The variational model sheds some light on this 
regime by clarifying the role of the temperature in systems out of equilibrium, but there is some way to go before we understand 
issues associated with, for example, any "principle of extremal entropy production" and instabilities that lead to structure for- 
mation. Finally, despite the obvious successes of the extended thermodynamics framework ilql . there is no universal agreement 
concerning the validity (and usefulness) of the results. To some extent this is natural given the interdisciplinary nature of the 
problem; to make progress we need to account for both thermodynamical principles and fundamental general relativity. This 
leads to a range of deep questions concerning, in particular, the actual meaning of the variables involved in the different models 
(e.g., the entropy). The ultimate theory must have a clear link with statistical physics and even information theory. Our efforts 
are not yet at that level. Basically, we need to continue to make progress if we are to address fundamental problems in, for 
example, cosmology. 

This paper sets a rather more modest target; we want to explore the extent to which a "first-order" formulation for heat 
conduction in general relativity is viable. The question may seem somewhat odd given that we have already acknowledged 
the need to account for (at least) second order contributions. However, it is interesting to ask whether a first-order model may 
nevertheless be useful (possibly in a somewhat restricted sense). We will demonstrate that this is, indeed, the case. Noting 
that the original first order models, due to Eckart [171 and Landau and Lifshitz [I18II . were incomplete we develop a consistent 



framework that includes the key thermal relaxation. We then consider the properties of this model, and show that it can be made 
both stable and causal (making contact with the classic work by Hiscock and Lindblom Ill9ll20ll ). This does not mean that the 
system may not exhibit instabilities. On the contrary, instabilities are in a sense generic in these problems 121!], but the analysis 
sheds further light on the nature of these instabilities and also elucidates the stabilizing role of the thermal inertia. The discussion 
also provides insight into the emergence of second sound (an effect that has been experimentally verified in low temperature 
crystals) associated with the heat transport. This provides a key link to systems that exhibit superfluidity, and demonstrates the 
potential for a unified treatment of heat transport in normal and superfluid matter. 

II. THERMAL DYNAMICS 

We take as our starting point our recent variational analysis of the relativistic heat problem [12] . The model is phenomenolog- 
ical, and assumes that the entropy component can be treated as a "fluid". In essence, this implies that the mean free path of the 
phonons is taken to be small compared to the model scale. We then consider two fluxes, one corresponding to the matter flow 
and one which is associated with the entropy. The latter is essentially treated as a massless (zero rest-mass) fluid. The dynamics 
then follows from a Lagrangian which depends on the relative flow of the two fluxes. The associated entrainment turns out to be 
a crucial feature of the problem 11121 [l4ll22|] . 

We assume that the particle number is large enough that the fluid approximation applies and there is a well defined matter 
current, n°. Moreover, we adopt the multi-fluid view and treat the entropy as an effective fluid with flux s". This current is in 
general not aligned with the particle flux. The misalignment is associated with the heat flux and leads to entropy production. 

As in the case of a generic two-fluid system (see |23] for an example in the case of a cool relativistic superfluid), the starting 
point is the definition of a relativistic invariant Lagrangian A. Assuming that the system is isotropic, we take A to be a function 
of the different scalars that can be formed by the two fluxes 113811 . From n" and s° we can form three scalars; 



n = —nnii 



J = -UaS 



(1) 



An unconstrained variation of A then leads to 



,, dA, dA, dA^. 
dA = -^On + -—ds + -—dj 
on OS O] 



(2) 



Changing the_passive density variations for dynamical variations of the worldlines generated by the fluxes and the metric (as 

dA 



discussed in f 1]) we find that 
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From this result we can read off the conjugate momentum associated with each of the fluxes; 

dn"- ds"- 



where we have introduced the coefficients IIU [T2ll 



B" = -2|4> ^' 



^2^ A" 
ds'' 



dA 



(4) 



(5) 



The energy-momentum tensor is obtained by noting that the displacements of the conserved currents induce a variation in the 
spacetime metric and therefore the variations of the fluxes, Sn" and ds", are constrained. The energy-momentum tensor is thus 
found to be 



Tj> = ^Ian'' + Oas" + ^S^' , 
where we have defined the generalized pressure, ^, as 

* = A - flan" - OaS" . 



(6) 



(7) 



As a result of the coordinate invariance associated with general relativity, the divergence of the energy-momentum tensor ^ 
vanishes. For an isolated system, we can express this requirement as an equation of force balance 



ybTj> = /r + /: - , 



(8) 



where the individual force densities are y 

f: = 2n'>V[,^la] + f-^a'^bu" , (9) 

/: = 2s'y[t0a] + Oa^bs' . (10) 

We note that, in order to obtain the energy momentum tensor ^ we needed to impose the conservation of the fluxes as 
constraints on the variation 1 1]. However, the equations of motion, ^ and ( fTol i. still allow for non-vanishing production terms. 
If we, for simplicity, consider a single particle species, the matter current is conserved and we have Va^" — 0. This removes 
the second term from the right-hand side of ^. In contrast, the entropy flux is generally not conserved. In accordance with the 
second law, we must have 

VaS"=r, >0. (11) 

A. Temperature 

To make progress, we need to connect the general variational results with the relevant thermodynamical concepts. In doing 
this it makes sense to consider a specific choice of frame. In the context of a single (conserved) species of matter, we see that 
force /" is orthogonal to the matter flux, n°, and therefore it has only three degrees of freedom. Furthermore, because of the 
force balance ([8]), we also have ?i°/^ = 0. This suggests that it is natural to focus on observers associated with the matter frame. 
We therefore introduce the four- velocity u° such that n°- = nu"-, where Uau"' = — 1 and n is the number density measured in 
this frame. 

Having chosen to work in the matter frame (in the spirit of Eckart [171]), we can decompose the entropy current and its 
conjugate momentum into parallel and orthogonal components. The entropy flux is then expressed as 

s" = s* {u" + w"-) , (12) 

where w"- is the relative velocity between the two fluid frames, and u°'Wa ~ 0. Letting s" ~ su^ where u^ is the four-velocity 
associated with the entropy flux, we see that s* = 57 where 7 is the redshift associated with the relative motion of the two 
frames 1 39] . This illustrates the subjective nature of entropy. It is an observer dependent quantity, not an absolute notion. 
Similarly, we can write the thermal momentum as 

da ^ 0*Ua + B''s*Wa = {B" S* + A^'fl) Ua + 6"s*W„ . (13) 

This leads to a measure of the temperature measured in the matter frame; 

- u''6l„ = 61* = S^s* + ^"'^n . (14) 

In essence, this quantity represents the effective mass of the entropy component. Returning to the stress-energy tensor, and 
making use of the projection orthogonal to the matter flux, we find that the heat flux (energy flow relative to the matter) is given 
by 

ga = - -Lab UcT'"' = S*9*Wa , (15) 

where we have used the projection 

±^b^ gC^b ^ U^'u'' . (16) 

Defining the new variables ct" — s*w'^ and pa — B^s*Wa, the energy density measured in the matter frame can be obtained by 
a Legendre-type transform on the master function. That is, we have 

p*-U„Ubr"''--A+PaCTV (17) 

The relevance of the new variables becomes apparent if we consider the fact that the dynamical temperature in (fl4] l agrees 
with the thermodynamical temperature that an observer moving with the matter would measure. In other words, we have 



dp* 
ds* 



(18) 



where p* = p*{n,s* ,p). This is, of course, the standard definition of temperature as energy per degree of freedom of the system. 
Mathematically, the temperature is obtained from the variation of the energy with respect to the entropy in the observer's frame 
(keeping the other thermodynamic variables fixed). 



This result is not trivial. The requirement that the two temperature measures agree determines the additional state parameter, 
p, to be held constant in the variation of p* . The importance of the chosen state variables is emphasized further if we note that, 
when the system is out of equilibrium, the energy depends on the heat flux (encoded in ct" and pa). This leads to an extended 
Gibbs relation (similar to that postulated in many approaches to extended thermodynamics jlql ): 

dp* = pdn + e*ds* +adp. (19) 

This result arises naturally from the variational analysis. 

According to the traditional view, thermodynamic properties like pressure and temperature are uniquely defined only in 
equilibrium. Intuitively this makes sense since, in order to carry out a measurement, the measuring device must have time to 
reach "equilibrium" with the system. A measurement is only meaningful as long as the timescale required to obtain a result is 
shorter than the evolution time for the system. However, this does not prevent a generalisation of the various thermodynamic 
concepts (as described above). The procedure may not be "unique", but one should at least require the generalised concepts to 
be internally consistent within the chosen extended thermodynamics model. Our model satisfies this criterion. 

B. Causal heat flow 

The variational model encodes the finite propagation speed for heat, as required by causality. To see this, we use the orthogo- 
nality of the entropy force density /° with the matter flux, solve for the entropy production rate Fg and finally impose the second 
law of thermodynamics. It is natural to express the result in terms of the heat flux g", defined by 

s" = s*u'' + ^q'^ . (20) 

We also let the conjugate momentum takes the form 

9a^0*Ua+/3qa, (21) 

where we have defined 

1 A^^'n" 



s*e* 



(22) 



With these definitions, we impose the second law of thermodynamics by demanding that the entropy production is a quadratic 
in the sources, i.e., 

r.^^>0, (23) 

where k > is the the thermal conductivity. This means that the heat flux is governed by 

T {(f + qcV^'u^) + g" = -K 1"'' {Vb9* + e*ub) , (24) 

where if = u^'Vbq'^ and u" is the four-acceleration (in the following, dots represent time derivatives in the matter frame). Here 
we have introduced 

(25) 



1 + K^ 



while the thermal relaxation time is given by 



kB 

(26) 



l + K/3 



The final result ( [24b is the relativistic version of the so-called Cattaneo equation 111211 . From the analysis we learn that the entropy 
entrainment, encoded in A'^^, plays a key role in determining the thermal relaxation time r. This agrees with the implications of 
extended thermodynamics, and echoes recent results in the context of Newtonian gravity Ill4ll . 



C. The matter flow 

The heat problem has two dynamical degrees of freedom. So far, we have focussed on the entropy. In addition to the relativistic 
Cattaneo equation ( l24b we have a momentum equation for the matter component. From ([9]) it follows that this equation can be 
written 

l-ma+ J-a Vh^J. + aqa + aqa + aq^VaUt = -/" . (27) 

n 

Here we have represented the matter momentum by 

fJ-a = ^J■Ua + aqa , (28) 

where /i is the chemical potential (in the matter frame) and 



This means that we have 



Given these definitions, we have [c.f., ( |23]) 1 



(29) 



a^'-^. (30) 



Ja ^ Ja ~ ~ 



s 



qa ■ (31) 



It is useful to note that this implies that the force has a term that is linear in g". This will be important later 
These relations complete our summary of the heat conduction model developed in II12I1 . 

III. A CONSISTENT FIRST-ORDER MODEL 

The model we have described crucially contains terms that enter at second order of deviation from thermal equilibrium, e.g., 
terms that are second order in the heat flux g". Moreover, it is clear that key effects (like the entropy entrainment) arise from 
the presence of second order terms in the Lagrangian A. Having said that, it is obviously the case that we can truncate the 
model at first order This does not take us back to the original first-order model discussed by Eckart [17.1 . Crucially, the thermal 
relaxation remains. Basically, this reflects the simple fact that you need to know the energy of a system to quadratic order in 
order to develop the complete linear equations of motion. Noting this, it is interesting to consider the features of this new first- 
order model. First of all, we can expect to get a clearer understanding of some of the general features of the variational model. 
Secondly, we may also find that this, much simpler, model is adequate for many situations of practical interest. 

A. The hnear model 

We want to restrict our analysis to first order deviations from equilibrium. Thermal equilibrium corresponds to g" = 0, no 
heat flux, and ?i" = 0, no matter acceleration. Moreover, in the simplest cases there should be no shear, divergence or vorticity 
associated with the flow, i.e., we will have Vqw" = and Vbu"" = as well |40]. Treating all these quantities as being of first 
order, and noting that 

Ubq'' = -q^iib , (32) 



contributes at second order, we arrive at two momentum equations; from (I27t we have 

^lUa+ ^''a '^bfJ- + aqa+ (a -) qa=0 , (33) 

while (|24ll leads to 

rqa +qa + fi {^\ ^bT + Tu^) = , (34) 



We also have the two conservation laws 

V^n'^ = , (35) 

Vas" = . (36) 

In these equations we have used the fact that s* and 6* differ from the equilibrium values s and T only at second order Moreover, 
to first order the pressure p is obtained from the standard equilibrium Gibbs relation; 

VaP = nVafi + sVaT . (37) 

Finally, we have the fundamental relation 

p + p = ^n + sT . (38) 

By comparing ( [33] l and ( |34] l to Eckart's results it becomes apparent to what extent the first-order model remains cognizant of its 
higher order origins. Specifically, a and (therefore) r depend on A'^^ and the entropy entrainment, c.f., ( |29] l. These effects rely 
on quadratic terms in the Lagrangian, and hence would not be present in a model that includes only first order terms from the 
outset. Hence, they are absent in Eckart's model. 

In order to analyze the dynamics of the heat problem, we will consider perturbations (represented by 6) away from a uniform 
equilibrium state. First of all, we have Qa — Ua ^ for a system in equilibrium. We can also ignore a and /?, since the 
equilibrium configuration is uniform, which means that we can replace k by k. This means that we are left with the two 
equations; 

^lSua+ 1-^ ^bS^i + aSqa - —Sqa = , (39) 

UK 



and 



rSqa + 6qa + K l.\ V^ST + nTSua = , (40) 



It is worth noting that we can combine these two to get 

{p + p) 5Ua+ l^a ^bSp + Sqa = . (41) 

The last two equations [ (|40| | and dTTT ll are, not surprisingly, identical to the first-order reduction of the Israel-Stewart model. This 
means that the problem is relatively well explored. In particular, the conditions required for stability and causality were derived 
by Hiscock and Lindblom 1 19, 20], see also Olson and Hiscock |25], quite some time ago. However, there are good reasons to 
revisit the problem. Most importantly, there is clear evidence from the recent literature (c.f., discussions of the relevance of the 
thermal relaxation and the role of the coupling between the four acceleration and the heat flux ll26l - l28ll ) that the key lessons from 
almost three decades ago have not been appreciated. To some extent this could be due to the fact that the Hiscock-Lindblom 
analysis is rather involved. Our aim is to clarify the main issues in the simpler context of heat conduction (ignoring viscosity). 
We also want to emphasize aspects that were only mentioned in passing in early work. Particularly relevant in this respect is 
the existence of second sound; an effect that is prominent in superfluids but which has also been observed in low-temperature 
crystals. We will demonstrate how the second sound emerges within the causal heat-conduction model. The overarching aim 
is to establish, beyond any reasonable doubt, that the model represented by ( |40] | and either ( l39b or (HTI) has all the properties 
expected of a reliable model for heat conduction in general relativity. 

B. Transverse waves: Stability 

Working in the frame associated with the background flow, we note that ( [39] ) and ( |40] | only have spatial components. That is, 
we may erect a local Cartesian coordinate system associated with the matter frame and simply replace a ^f i where i ^ 1 — "i. 
Then taking the curl (e^'^* Vfc) of the equations in the usual way, we arrive at 

m.U' - ig* = , (42) 

r 

and 

TO*Q' + {p + p)Q' = , (43) 



where we have defined 

U' = e'l'^V j5uk , and Q' = e'^^V j5qk , (44) 

and 

/ aKT\ kT 
rrii, = nlfi 1 = p + p . (45) 

Assuming that the perturbations depend on time as e'"*, where t is the time-coordinate associated with the matter frame, we 
arrive at the dispersion relation for transverse perturbations; 

ioj[{p + p){l+iojT) -iojKT]=Q . (46) 

Obviously w = is a solution. The second root is 

i(p + p) 
uj = -Si--J2. . (47) 

This results shows that the thermal relaxation time r is essential in order for the system to be stable. We need m* > 0, i.e., the 
relaxation time must be such that 

kT 
T > —— . (48) 

p + p 

The analysis clearly shows that Eckart's model (for which t = 0) is inherently unstable. Moreover, the constraint on the 
relaxation time agrees with one of the conditions obtained by Olson and Hiscock |251 (c.f., their eq. (41)), representing the 
inviscid limit of the exhaustive analysis of the Israel-Stewart model of Hiscock and Lindblom 111911 . We may also note that the 
condition given in eq. (43) of |25] simply leads to the weaker requirement r > 0. 

The physical origin of the transverse instability can be understood if rewrite ( l40l i and ( l4Tl i as 

mjua+ ±i VbSp - i {Sqa + K ±!i VbST) = , (49) 

and 

mjqa + {p + p) {Sqa + k l''^ ST) - kT l^ VbSp = . (50) 

These relations show that to* plays the role of an "effective" inertial mass (density). The importance of this quantity has been 
discussed in a series of papers by Herrera and collaborators II29I - DII1 . especially in the context of gravitational collapse. Basically, 
the instability of the Eckart formulation is due to the inertial mass of the fluid becoming negative. Once this happens the pressure 
gradient no longer provide a restoring force, rather it would tend to push the system further away from equilibrium. This is a 
run-away process, associated with exponential growth of the perturbations. Ultimately, the instability is due to the inertia of 
heat; an unavoidable consequence of the equivalence principle (heat carries energy, which means that it can be associated with 
an effective mass 11321] '). The condition ( l48T l may seem rather extreme (Hiscock and Lindblom |20] quote a timescale of 10^^^ s 
for water at 300K), but it sets a sharp lower limit for the thermal relaxation in physical systems. A system with faster thermal 
relaxation can not settle down to equilibrium. However, it may still be reasonable to ask if a system may evolve in such a way 
that it enters the unstable regime (in the way discussed in ll29ll30ll ). Given our assumed equilibrium the present formulation does 
not allow us to consider this question, but it seems clear that if a system were to evolve in that way then one would need a fully 
nonlinear analysis to determine the consequences. 

C. The longitudinal problem 

The transverse problem is relatively simple since the are no corresponding restoring forces in a simple fluid problem (these 
requires rotation, elasticity, the presence of a magnetic field etcetera). When we turn to the longitudinal case the situation 
changes. In a perfect fluid longitudinal perturbations propagate as sound waves, and when we add complexity to the model the 
dispersion relation can get very complicated. 

To analyze the longitudinal case, we take the divergence of (|40] | and PTt . This leads to 

{p + p)dtVM' + y^Sp + dtVM' = , (51) 



and 

TdtViSq' + ViSq' + kV^ST + KTdtV^Su' = . (52) 

We also need the conservation laws which take the form 

VM' ^--dtSn, (53) 

n 

ViSq' = -nTdtSs , (54) 

where we have defined the specific entropy; s = s/n. To make progress we need to, first of all, decide what variables to work 
with, and secondly we need an equation of state for matter In the following we will opt to work with the perturbed densities 6n 
and Ss. Keeping in mind that we are only retaining first-order quantities, we have 

*^(I)/-(S),/- 

and 

It is also useful to keep in mind that the temperature represents the "entropy chemical potential", i.e. is defined by 

where p = p{n, s) represents the equation of state. This immediately implies that 

dn J ^ \ ds 



(58) 



i.e., we can reduce the number of thermodynamic quantities. In order to make contact with (potential) observations, it is natural 
to work with (i) the adiabatic speed of sound; 



2 i ^P\ "- f ^P 



(ii) the heat capacity at fixed volume; 



and (iii) 



''^-[fp)rjT-p[tn) ' ^^^) 



c. = r(^) , (60) 



n fdT\ T (dp\ 



For future reference, it is also useful to note the identity [c.f., eq. (96) in Hiscock and Lindblom |1 

1 1 _ n^ fdTV _ nT 2 

Cy Cp T{p + p)c2 \dnjg {p + p)cl " ' 

where Cp is the heat capacity at fixed pressure. 

If we (again) focus on plane-wave solutions such that the perturbations behave as cxp(ia;i + ikjX^ ), and introduce the phase- 
velocity a = w/fc, then the above relations lead to the coupled equations 

^-i-^ (cr^ - cl) 6n + nT (a^ - as) Ss = , (63) 



and 



(tr — as) on + 1 nra 

n V k Cr, 



5s ^Q 



(64) 



From these results it is easy to see that the generic dispersion relation will be a quartic in a. This is exactly what one would expect 
for a physical system. As we have already mentioned, the formulation must accommodate the presence of "second sound" which 
has been experimentally verified for both superfluids and low temperature solids. Working out the explicit dispersion relation, 
we find 



iijp + p) 



^{-'-cD- 






c,.T — 2KTa. 



P + PCs 



Tat 



n 



= 



(65) 



This expression is still too complicated for us to be able to make definite statements about the solutions, without making further 
assumptions. The most direct strategy would be to consider an explicit equation of state, work out the relevant thermodynamics 
quantities, and then solve ( l65T l (probably numerically). This would allow us to establish whether the considered model is stable 
and causal. This route is, however, not particularly attractive given the need to introduce an explicit model. If we want to 
continue to consider the problem in (at least to some extent) generality, then we need to resort to approximations. As we will 
see, this is a very instructive route. 

In order to simplify the analysis, we will consider the long- and short-wavelength limits of the problem. The results we obtain 
in these limits provide useful illustrations of the key features. At the same time, we should keep in mind that both cases are 
somewhat "artificial". First of all, fluid dynamics is, fundamentally, an effective long-wavelength theory in the sense that it arises 
from an averaging over a large number of individual particles (constituting each fluid element). In effect, the model only applies 
to phenomena on scales much larger than (say) the interparticle distance. However, the infinite wavelength limit represents a 
uniform system, which is artificial since real physical systems tend to be finite. Moreover, as we will not account explicitly for 
gravity we can only consider scales on which spacetime can be considered flat. While the plane-wave analysis holds on arbitrary 
scales in special relativity, a curved spacetime introduces a cut-off lengthscale beyond which the analysis is not valid. The 
theoretical framework is valid in general, but on larger scales we would have to consider also the perturbed Einstein equations. 



D. The long wavelength problem 



Let us first consider the long wavelength, fc ^ 0, problem. This represents the true hydrodynamic limit, and it easy to see that 
there are two sound-wave solutions and two modes that are predominantly diffusive. The sound-wave solutions take the form 



±Cs 



l±i 



kT 



2(p + p)c3 



{cl-as^k 



These solutions are clearly stable, since Im cr > 0. Using the Maxwell relations listed by Hiscock and Lindblom y. 
show that this results agrees with eq. (40) from ll20ll . Moreover, our result simplifies to [using ( |62] l1 



(66) 



we can 



Im a 



K 

2n 



1 



1 



(67) 



in the limit where \as\ ^ c^, which is relevant since c^ ^ p/p becomes small in the non-relativistic limit. Indeed, we find that 
(l67T l agrees with the standard result for sound absorption in a heat-conducting medium IBSll . 
In addition to the sound waves, we have a slowly damped solution 



1 



Tai 



This is the classic result for thermal diffusion. 
Finally, the system has a fast decaying solution; 



{p + p)c; 



i{p + p) 

m^kr 



(68) 



(69) 



Under most circumstances, this root decays too fast to be observable. This means that the model reproduces that standard 
"Rayleigh-Brillouin spectrum" with two sound peaks symmetrically placed with respect to the broad diffusion peak at zero 
frequency |;33;,i34;l 
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E. Short wavelength stabihty and causality 



Different aspects of the problem are probed in the short wavelength limit. Letting /c — > oo we see 114 111 that (1651 ) reduces to a 
quadratic for a^ . This allows us to write down the solutions in closed form. Hence, it is relatively straightforward to establish 
the conditions required for the stability of the system in this limit. For infinitesimal wavelengths we have 



1 



2A 



Acr'^-Ba^ +C = 0, — > cri 

where 

A ~ m^.T > 
in order for transverse perturbations to be stable. We also have 



B±{B^ - 4:AC) 



1/2- 



B = {p + p) 



cIt 



2kTq;, 



and 



C = K 



p + pel 



Tai 



^>0. 

Cp 



(70) 



(71) 



(72) 



(73) 



In order to guarantee stability for longitudinal perturbations, we need a^ to be real and positive. Given the quadratic formula 
and the fact that A> Q this implies that we must have B^ — AAC > 0. After some algebra, this leads to 



C^.T 



K 



2KTas 
P + P 



P + P 



kT 



Ak'T 



[p + p)ncv 



- 2a,) > . 



(74) 



The first two terms are positive, as long as ( |48] | is satisfied. Hence, the condition is guaranteed to be satisfied as long as c^ > as- 
As discussed later, this is certainly the case for degenerate matter. In cases where this simple condition is not satisfied, ( |74] l 
provides a complicated constraint on the relaxation time. Finally, we must also have i? > 0, which can be expressed as; 



T > 



2r 

p + p 



-a. 



1 

71C„ 



(75) 



and l//i?o and l//?2 both also 



This condition is identical to that given in eq. (146) of 111911 (obtained in the limit where Ui 
vanish, c.f. jHHl) 

Let us now consider finite wavelengths. Letting a ~ (7±+ <Ji/k, where a± solve (fTOl l. and linearising in 1/k, we find that 



CTl 



i{p + p) ( o\-cl 



2 V 2^cr^ - B 



(76) 



Since all quantities in this expression are already constrained to be real, we need Im ci > (for real fc) in order for the system 
to be stable.. From (|70] i we see that 



2Aa\ 



This then leads to the final condition 



B = ±LB^ -4AC 



2^2^ 2 

cr_ < c„ < (T I 



1/2 



(77) 



(78) 



It is worth noting that this result is consistent with the notion that "mode-mergers" signal the onset of instability f21|l. 

As the waves in the system must remain causal, we should insist that a^ < 1. To ensure that this is the case, we adapt the 
strategy used by Hiscock and Lindblom ll9ll . As JTOl i is a quadratic for a^ we can ensure that the roots are confined to the 
interval < cr-^ < 1 (noting first of all that the roots are real since (|74] | is satisfied). Given the B and C are both positive, the 
roots must be such that a^ > 0. Meanwhile we can constrain the roots to cr^ < 1 by insisting that 



A-B + C>0, 



(79) 



and 



A-2B>0. 



(80) 
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Combining these inequalities with the positive discriminant, we can show that A > B/2 > C. The first of the two conditions 
can be written 



(1-a 



K 

T 

ncv 



2 



^«T(i^a^^O (81) 

P + P 



Now, when combined with causaUty the condition (|78] l requires that c^^ < (j\ < 1. In other words, we must have c^ < 1, which 
means that dSTT l implies that 

T > — . (82) 

Comparing to the results of Hiscock and Lindblom lll9ll . we recognize (IHTj ) as their fia > condition (it is also eq. (4) of Herrera 
and Martinez IBOtl ). while ( l82b corresponds to fig > 0. 
Meanwhile, the condition dSOl l can be written 

K 2kT 

{2~cI)t> + (1-a,), (83) 

nc.„ p + p 

corresponding to eq. (148) of Hiscock and Lindblom. Finally, A > C leads to 

r > \ . (84) 

p + p ncp 



This corresponds to eq. (3) in Herrera and Martinez [30], which derives from eq. (147) of Hiscock and Lindblom JlS 

This completes the analysis of the stability and causality of the system. We have arrived at a set of conditions on the thermal 
relaxation time (and related them to results in the existing literature). As long as these conditions are satisfied, the solutions to 
the problem should be well behaved. 

F. The emergence of second sound 

So far we have considered the conditions that must be satisfied by a thermodynamical model in order to ensure the stability 
and causality of both transverse and longitudinal waves. To complete the analysis of the problem, we will now consider the 
nature of the solutions. Since the phase velocity a is obtained from a quartic, we know that the problem has two (wave) degrees 
of freedom. This accords with the experience from superfluid systems and experimental evidence for heat propagating as waves 
in low temperature solids. One of the solutions should be associated with the usual "acoustic" sound while the second degree of 
freedom will lead to a "second sound" for heat. We want to demonstrate how these features emerge within our model. 

In order to explore these features, it is natural to consider the large relaxation time limit. Taking the relaxation time r to be 
long, the solutions to ( fTOl i take the form (up to, and including, order 1/r terms) 

(85) 
L (P + P}T \ ciJ \ 

which could be rewritten using (l62b . and and 

(86) 

riTCp 

The first of these solutions clearly represents the usual sound, while the other solution provides the second sound. In the latter 
case, the deduced speed is exactly what one would expect 116:1 . It is easy to see that the first root will satisfy jTSl ), and the 
associated roots will be unstable in the long relaxation time limit. Moreover, the second solution leads to stable roots provided 

T > ^^ . (87) 

Basically, the finite wavelength condition implies that the second sound must propagate slower than the first sound. This is, 
indeed, what is measured in physical systems (like superfluid Helium). Moreover, it is easy to see that this condition must be 
satisfied in order for the long relaxation time approximation to be valid. 

As a useful illustration of the properties of the model, let us consider the particular case of degenerate matter In this case, 
which relates to electrons in both metals and white dwarfs (and also neutrons and protons in neutron stars), the two specific heats 
are almost identical; 

^«l + ofM:y, (88) 

Cp \ ^F ) 



\ 


2 


kT 

1 + 7 ^ 

[P + P)r 


(-1)] 
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FIG. 1 : This figure provides an illustration of the qualitative nature of the behaviour of heat conducting degenerate matter (as discussed in the 
main text), based on the consistent first-order relativistic model. The parameters have been chosen in such a way that the speed of sound is 
10% of the speed of light, while the second sound (at short wavelengths, large k) propagate at l/\/3 of this. The phase velocity of the waves is 
a = Re Ld/k (left panel). The thermal relaxation time r has been chosen such that the critical wavenumber at which the second sound emerges 
is fc = 10. At lengthscales larger than this, the corresponding roots are diffusive (have purely imaginary frequency), and in the very long 
wavelength limit (A: — ^ 0) we retain the expected thermal diffusion. The damping time follows from 1/Im uj (right panel).We also indicate the 
acausal region (grey are). The illustrated example is clearly both stable and causal. 

where ks is Boltzmann's constant and ep is the Fermi energi [36]. This means that, for temperatares significantly below the 
Fermi temperature, we can accurately assume that as ~ 0. If we also assume that t ^ nT/{p + p) then the dispersion relation 
factorises and we have 



1 



K 

ncv 



(89) 



That is, the four roots are 



a = ±Cs 



(90) 



and 



2kT 



l±(l-^fc 



1/2- 



(91) 



The character of these roots is illustrated in Figure [T] We see that the ordinary sound exists at all wavelengths. Meanwhile, at 
short long wavelengths (small fc) the remaining two roots are exponentially damped, i.e. diffusive in character. One root has a 
relatively slow decay, corresponding to the expected thermal diffusion, while the other root decays so rapidly that it is unlikely 
to be observable by experiment. Below a critical lengthscale (corresponding to A: = 10 in Figure [T]) the second sound emerges 
as a result of the finite thermal relaxation time r. For very short lengthscales, heat signals will propagate as waves. However, as 
is clear from ( |9T1 ) these solutions are always damped. In order to 'propagate' the real part of the wave frequency must exceed 
the imaginary part (so that several cycles are executed before the motion is damped out). This boils down to the second sound 
propagating only for 



nCv 



-1/2 



(92) 



This result is interesting if we consider systems that become superfluid (see [37] for an interesting discussion of models for 
relativistic superfluids). Suppose we consider a system which starts out in the diffusive regime (e.g. Helium above the superfluid 
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transition temperature). When the system is cooled down through the relevant transition temperature, the (non-momentum 
conserving) particle collisions that give rise to k are suppressed. In effect, the critical value of k decreases and the system may 
enter the regime where the second sound can propagate at macroscopic scales. The change to the basic nature of heat propagation 
is easily understood. 

At the end of the day, this relatively simple model demonstrates how the behaviour of a given physical system depends on the 
balance between the characteristic timescales. Obviously, we also need to keep in mind that real systems impose restrictions both 
for large k, as the fluid model breaks down when we approach the interparticle distance scale, and small k, when the wavelength 
becomes larger than the size of the system. 

IV. FINAL REMARKS 

We set out to derive a consistent first-order model for relativistic heat conduction, in such a way that the theory remains 
cognizant of its higher order origins. As we have demonstrated, this leads to a model that retains the thermal relaxation that is 
necessary if we want the problem to remain causal. What have we learned from this exercise? First of all, we have illustrated 
the problems associated with the original first-order models, e.g. that of Eckart 1.17.1 . The conclusions are, obviously, not new 
but the discussion should lay to rest any suggestions that the coupling of the four-acceleration to the heat-flux is (somehow) 
problematic [26-28]. By considering the waves present in the new model, we have established that the system is both stable and 
causal provided that some seemingly natural conditions are satisfied. This conclusion accords with the classic work by Hiscock 
and Lindblom Ill9ll20ll (see also Olson and Hiscock ||25I1 ). In fact, the conditions we have arrived at reproduce their key results. 
However, our analysis adds to previous work by discussing the emergence of the second sound and the nature of the associated 
solutions. This is a key point, especially if we are interested is relativistic superfluid systems. The analysis also demonstrates the 
intricate nature of these problems. It is easy to see how a model that may fail one, or several, of the derived conditions in some 
regime may nevertheless be valid for a different range of parameters. Hence, one really should consider the applicability of the 
chosen theory on a case-by-case basis. This is probably no more than should be expected from a phenomenological model. 

Our results represent useful progress in this problem area, but one could obviously develop the theory further A natural step 
would be to consider the various constraints that we have derived for detailed equations of state, e.g., matter coupled to phonons. 
It would also be interesting to consider applications of the first-order construction. While the model is restricted in the sense that 
it does not account for non-adiabatic effects, there is an exciting range of possible applications in astrophysics, cosmology and 
high-energy physics. Particularly interesting questions concern to what extent second sound effects are relevant in relativistic 
systems and the difference between first-order results and the, considerably more complex, second-order theories. 
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